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The H.E.S.S. array of Cherenkov telescopes has performed, from 2004 to 2007, a survey of the inner 
Galactic plane at photon energies above 100 GeV. About 400 hours of data have been accumulated 
in the region between -30 and +60 degrees in Galactic longitude, and between -3 and +3 degrees in 
Galactic latitude. Assuming that dark matter is composed of Weakly Interacting Massive Particles, 
we calculate here the H.E.S.S. sensitivity map for dark matter annihilations, and derive the first 
experimental constraints on the "mini-spikes" scenario, in which a gamma-ray signal arises from 
dark matter annihilation around Intermediate Mass Black Holes. The data exclude the proposed 
scenario at a 90% confidence level for dark matter particles with velocity-weighted annihilation 
cross section av above 10 -28 cm 3 s _1 and mass between 800 GeV and 10 TeV. 

PACS numbers: 98.70.Rz, 98.56.Wm, 95.35.+d 
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I. INTRODUCTION 

A substantial body of cosmological and astrophysi- 
cal measurements suggests that ^22% of the Universe 
is composed of non-baryonic dark matter (DM), e.g. 
commonly assumed to be in the form of Weakly Inter- 
acting Massive Particles (WIMPs) arising in extensions 
of the Standard Model of Particle Physics (for recent re- 
views see @, Q ) ■ The lightest neutralino arising in super- 
symmetric (SUSY) extensions of the Standard Model [4j, 
and the first excitation of the Kaluza-Klein bosons,!^ 1 ), 
in universal extra dimension (UED) theories 0,0, 0] are 
among the most widely discussed dark matter candi- 
dates. 

Although accelerator and direct dark matter searches 
may well provide useful hints on the nature of dark 
matter, a correct identification is likely to require the 
combination of different techniques, including the indi- 
rect searches, based on the detection of the annihila- 
tion products of dark matter particles. The annihila- 
tion rate being proportional to the square of the dark 
matter density integrated along the line of sight, re- 
gions with enhanced dark matter density are primary 
targets of indirect searches. Among them are the Galac- 
tic halo 0, external galaxies p i, gala xy clusters [l(|, 
substructures [ill HI [lE flj [TjJiaMI [H ll^ ■ and the 
Galactic Center (GC) 0, E|, El, [3 . 

The GC, in particular, has attracted significant inter- 
est. The distribution of dark matter at the GC is ac- 
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tually highly uncertain, due to lack of resolution in N- 
body simulations, and to the many astrophysical effects 
that further complicate the situation, such as the pres- 
ence of the supermassive black hole coincident with Sgr 
A*, gravitational scattering of dark matter off the stellar 
cusp, and dark matter annihilation [25j |. 

Dwarf spheroidal galaxies in the Local Group (26| have 
also been considered as targets for gamma-ray detection 
since they seem to represent dark matter dominated re- 
gions [11 [H, H El [Hi [32|. More recently, the Sloan 
Digital Sky Survey (SPSS) revealed the existence of new 
satellites [H, [HI HI] offering appealing features for dark 
matter searches. Radial velocity dispersion of stars in 
Galaxy satellites usually implies large mass-to-luminosity 
ratio. Nevertheless, the lack of accurate measurements 
on the velocity dispersion may induce large systematic 
effects on the parameters used for the modelling of the 
dark matter halo. Even in case of accurate kinematic 
data for some of the Galaxy satellites, only faint annihi- 
lation signals are expected for smooth dark matter halos 
due to their distance of ~100 kpc. 

Mini-spikes around Intermediate Mass Black Holes 
(IMBHs) have been recently proposed as promising tar- 
gets for indirect dark matter detection [36}. Accurate 
predictions, in the context of well-defined astrophysi- 
cal scenarios, have been derived for the distribution and 
luminosity of these objects [37|- Mini-spikes might in 
fact be detected as bright point-like sources in gamma- 
rays [H| and neutrinos (39|, and the prospects for de- 
tection with satellites with lar ge fi eld of view such as the 
upcoming GLAST experiment [4tj] to be launched in 2008 
appear particularly promising. 

Current Imaging Atmospheric Cherenkov Telescopes 
(IACTs) can also effectively search for these objects. The 
H.E.S.S. (High Energy Stereoscopic System) experiment 
has already surveyed a significant part of the Galactic 
plane. With the combination of a large field of view, 
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very good angular resolution and off-axis performance, 
H.E.S.S. has reached the sensitivity to accurately map 
the Galactic plane in scan-based observations. 

In this paper, H.E.S.S. data are used to derive for the 
first time experimental exclusion limits on the dark mat- 
ter annihilation signals within the context of the mini- 
spike scenario. The paper is organized as follows. Sec- 
tion II is devoted to the mini-spike scenario and the 
gamma-ray flux expected from dark matter annihilations 
in mini-spikes. In section III, we present the H.E.S.S. 
data from the Galactic plane survey and compute the 
H.E.S.S. flux sensitivity map to dark matter annihila- 
tions in the Galactic plane region. Exclusion limits are 
then derived for dark matter annihilation from mini- 
spikes. Finally, section IV is devoted to the discussion 
of the results obtained in this study. 



II. DARK MATTER ANNIHILATIONS IN 
MINI-SPIKES 




A. Dark matter candidates in MSSM and 
Kaluza-Klein models 

In this paper we focus on two particle physics sce- 
narios beyond the Standard Model which provide well- 
motivated WIMP dark matter candidates with masses 
and couplings at the electroweak scale to account for the 
non-baryonic dark matter. The annihilation of WIMP 
pairs can produce in the final state a continuum of 
gamma-rays whose flux extends up to the DM particle 
mass, from the hadronization and decay of the cascading 
annihilation products. 

Minimal Supersymmetric extensions of the Standard 
Model (MSSM) predict the lightest supersymmetric par- 
ticle (LSPVto be stable in the case of R-parity conserving 
scenarios [J]. In various SUSY breaking schemes, this 
fermionic particle is the lightest neutralino \- I n gen- 
eral MSSM, the gamma-ray spectrum from neutralino 
annihilation is not uniquely determined and the branch- 
ing ratios (BRs) of the open annihilation channels are 
not determined since the DM particle field content is 
not known. Soft and hard spectra for the neutralino 
pair annihilation are therefore considered. The annihi- 
lation spectra are parametrized using PYTHIA [4l| sim- 
ulations. Fig. [TJ shows the number of gammas above the 
energy threshold E t h= 100 GeV, N 1 (E 1 > Eth), as a 
function of the neutralino mass m^. N 1 {E~ I > £th) is 
lower in the case of a 100% BR bb channel than in the 
t + t~ channel for m x below 2 TeV. Above 2 TeV, more 
gammas are expected in the bb channel. 

In some specific scenarios, the branching ratios of the 
annihilation channels can be computed given that the 
field content of the DM particle is known. In the SUSY 
scenario dubbed AMSB [H 0, H EH (Anomaly Me- 
diated Supersymmetry Breaking) arising in SUSY theo- 
ries 0, [13] , the DM candidate is the lightest neutralino 
with a predominant wino component which annihilates 



Figure 1: Number of gammas per annihilation above the en- 
ergy threshold E t h = 100 GeV, N 1 {E 1 > -Eth), as a function 
of the neutralino mass, m x , for 100% BR in the bb channel 
(solid black line) and 100% BR in the t + t~ (solid red line). 
This number corresponds to the annihilation spectrum inte- 
grated from E t h up to m x . Below 2 TeV neutralino masses, 
more gammas per annihilation are expected in the case of 
the t + t~ channel than for bb, thus the latter will provide a 
conservative value. Above 2 TeV, the bb channel yields more 
gammas than the t + t~ one. 



in W + W~ pairs decaying in quarks and leptons. Large 
annihilation cross sections are ex pec ted which may lead 
to detectable gamma-ray fluxes [4g. Within this sce- 
nario, neutralino masses may extend up to tens of TeV. 
The other dark matter candidate considered here arises 
in theories with universal extra dimensions (UED). In 
Kaluza-Klein (KK) scenarios with KK-parity conserva- 
tion, the lightest KK particle (LKP) is stable. Most of- 
ten, the LKP is the first KK mode of the hypercharge 
gauge boson, referred hereafter to as B^ 0|- In 

the Kaluza-Klein case, the BRs are extracted from Q. 

pairs annihilate mainly into fermion pairs: 35% in 
quark pairs and 59% in charged lepton pairs. 



B. IMBHs formation scenarios 

IMBHs are compact objects with masses comprised be- 
tween that of the heaviest remnant of a stellar collapse, 
~ 20M Q [i^, H3| , and the lower end of the mass range of 
Supermassive Black Holes (SMBH) - 1O 6 M (UHl. 

Observational hints of the existence of IMBHs come 
from the detection of ultra-luminous X-ray sources 
(ULXs), apparently not associated with active galactic 
nuclei (53l . l5~il [55| . from stellar kinematics in globular 



4 



-?.10 4 


A 



X 

CQ 



10 



3 =! 



10 



10 - 



2 _ 



1 1 1 1 1 1 



TTTTT] 1 | 1 1 ||U 



:>■„. 



crv = 3x10" 26 cm 3 s- 1 
— i — rrij = 1 TeV 
..... 111., = 500 GeV 
m, = 300 GeV 



i 

i i 



'A 

'. Ml 



10 18 10 17 10 16 10 15 10 14 10 13 10 12 „ 10 11 
<!> (E >100 GeV) [cmV] 





A 



X 

CQ 



102 f mmm H!!i|||jiL 



10 r 



1 1 1 1 1 1 



cv = 3x10~ 26 cmV 
—i — m z = 1 TeV 
.... 111, = 500 GeV 
m. = 300 GeV 




jQ-15 10 -14 jQ-13 1Q .12 10 11 jQ-10 10 -9 10 8 

<I> Y (E >100 GeV) [cm"Vl 



Figure 2: Integrated luminosity function of IMBHs for the scenario A, i.e. for IMBH masses of ~ 10 2 M (left), and the scenario 
B, i.e. for IMBH masses of ~ 10° M© (right), for a Milky Way-sized halo obtained from an average of 200 stochastic realizations 
(see text for details). Neutralino masses of 300, 500 GeV and 1 TeV, respectively, and a velocity- weighted annihilation cross 
section av = 3 x 10~ 26 cm 3 s _1 are considered. The integrated luminosity corresponds to the number of black holes yielding 
a integrated gamma-ray flux higher than a given integrated flux, Nbh(> $7), as a function of the integrated flux <E> 7 . The 
gamma-ray flux is integrated above 100 GeV. The nominal 5a H.E.S.S. point source sensitivity (25 hours) is plotted (black 
dashed-dotted line) for comparison. 



clusters (561 l57| and emission-line time lags in galax- 
ies [11] . 

From a theoretical point of view, a population of mas- 
sive seed black holes could help to explain the origin 
of SMBHs, especially those powering quasars at red- 
shift 6 or 7, that were thus already in place when the 
Universe was only ~ 1 Gyr old [591 ]. A population of 
"wandering" IMBHs is in fact a generic prediction of sce- 
narios that seek to explain the properties of supermassive 
black holes [6(J [Ml H2| • However, despite their theoret- 
ical interest, it is difficult to obtain conclusive evidence 
for the existence of IMBHs. 

In Ref. (38| the consequences of the formation and 
growth of IMBHs on the surrounding distribution of DM 
have been studied. In particular, it was shown that these 
processes lead to the formation of strong dark matter 
overdensities called mini-spikes, which are ideal targets 
for indirect dark matter searches, as they would appear 
as a population of gamma-ray point-sources with identi- 
cal energy spectrum. 

The properties of mini-spikes have been discussed in 
detail for two different scenarios. In the first one (sce- 
nario A), black holes are remnants of the collapse of 
Population III (or "first") stars [63], which are believed 
to collapse directly to black holes in the mass range 
M ~ 60 - 140 M Q and M > 260M Q (64|. Black holes in 
this scenario may not necessarily form at the very centers 
of their initial host dark matter halos at high redshift, a 



circumstance that, as we shall see, may have important 
consequences on the detectability of IMBHs. 

The second scenario (scenario B) is representative 
of a class of models in which black holes originate 
from the collapse of primordial gas in early-forming ha- 
los m m, m, m m, M, [at Thc initial black holcs 

are massive (~ 10 5 M Q ) and the growth of SMBH pro- 
ceeds in such a way that both mergers and accretion play 
an important role. Following Ref. [3o |. we focus here 
on the specific formation scenario proposed in Ref. [62| . 
which goes as follows: during the virialization and col- 
lapse of the first halos, gas cools, collapses, and forms 
pressure-supported disks at the centers of halos that are 
sufficiently massive to contain a relatively large amount 
of molecular hydrogen. Local gravitational instabilities 
in the disk lead to an effective viscosity that transfers 
mass inward and angular momentum outward [7l| until 
supernovas in the first generation of stars heat the disk 
and terminate this process [r32l ] . By the time the process 
terminates (of the order of the lifetimes of Pop III stars, 
~ 1 — 10 Myr), a baryonic mass of order ~ 10 5 A1© has 
lost its angular momentum and has been transferred to 
the center of the halo. Such an object may be briefly 
pressure-supported, but it eventually collapses to form a 
black hole (gllzi. 

The characteristic mass of the black hole forming in a 
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halo of virial mass M v is given by [6/ 



M bh = 3.8 x 1O 4 M 



3/2 



M„ 



1O 7 M 



■1+z 
, 18 



0.5) \0.Q3 J 

,3/2 / « \ 
) \ lOMyr J ' 



(1) 



where / is the fraction of the total baryonic mass in the 
halo that falls into the disk, z is the redshift of formation, 
K is the fraction of the baryonic mass which loses its 
angular momentum and remains in the remnant black 
hole, and t is the timescale for the evolution of the first 
generation of stars • 

In order to study the consequences for indirect dark 
matter searches we will compare H.E.S.S. data with the 
mock catalogs of Ref. [38| , which consisted of 200 stochas- 
tic realizations of Milky Way-like halos at z — 0, obtained 
by populating halos with black holes at high redshift (fol- 
lowing the prescriptions of scenarios A and B) and evolv- 
ing them forward to determine the properties of satellite 
black holes. The mass of the Milky Way was fixed at 
lO 12 ' 1 /?. -1 Mq at z = 0, and our analysis is based on 
their statistically large sample of wandering black hole 
populations in Milky Way-like halos of this mass. 

In scenario A, the mass spectrum of unmerged black 
holes is a delta function and the average number of un- 
merged black holes per Milky Way halo is Nbh,A — 1027± 
84, where the error bar denotes the la scatter from halo- 
to-halo [99]. In scenario B, the total number of unmerged 
black holes per Milky Way halo is iVbh.B - 101 ±22 
The dispersion denotes the la scatter among different 
realizations of Milky Way-like halos, as discussed in (38[. 



C. Mini-spikes 

The growth of massive black holes inevitably affects 
the surrounding distribution of Dark Matter. The profile 
of the final DM overdensity, called mini-spike, depends on 
the initial distribution of matter, but also on astrophysi- 
cal processes such as gravitational scattering off stars and 
mergers. 

Ignoring astrophysical effects, and assuming adiabatic 
growth of the black hole (i.e. assuming that the black 
hole grows on a time scale much longer than the dy- 
namical time scales of DM around it), one can calcu- 
late analytically the functional form of the final DM pro- 
file. If one starts from an initially uniform DM distri- 
bution, which is the most likely situation for black holes 
in scenario A, the final profile will be a mild mini-spike 
with density p sp oc (r/r/j) 3 / 2 (e.g. see [73| and refer- 
ences therein) . If one starts from a cuspy profile that is 
a power-law with index 7=1, as relevant for scenario B, 
the new profile is a new power-law, 



Ps P (r) = p(r sp ) — 



Sp 



(2) 



is related to the initial power-law index 7 by \7Z 

9-2 7 



7s 



7 



(3) 



In the case of the Navarro, Frenk and White profile, 7 = 
1, which implies 7 sp = 7/3. 

To calculate the annihilation flux, the singularity of 
p sp at r=0 needs to be cut off; we introduce a minimal 
radius r cut . One limit is given by the size of the IMBH, 
another by the condition that the annihilation rate of the 
dark matter particles is smaller than the inverse age of 
the mini-spike : 

Psp(nim) = m x /av(t-t f ) = p\ im . (4) 

An inner cut-off is therefore defined at a radius 



Max[4_R S chw,nim] 



(5) 



where i?schw is the Schwarzschild radius of the IMBH 
-Rschw = 2.95 kmAfbh/M©. For common values of the 
mass and cross section of the DM particle, r\i m ~ 10~ 3 pc 
so that r cut = n im . 

Although in principle this applies also to the black hole 
at the Galactic Center, there are a number of astrophys- 
ical effects, such as off-center formation, major mergers, 
and gravitational scattering off stars, that tend to erase 
any DM overdensity. A detailed discussion of the for- 
mation and evolution of the DM spike at the Galactic 
center, including a discussion of the prospects for indi- 
rect detection in light of the very high energy (VHE, 
E 7 > 100 GeV) gamma ray source coincident with Sgr 
A*, can be found in Ref. [251 ]. 

All these astrophysical processes are unlikely to take 
place around IMBHs. Mini-spikes around IMBHs that 
never experience mergers are therefore expected to be 
stable structures over cosmological timescales, and they 
are thus promising targets for indirect detection. The 
gamma-ray flux from these targets can be easily calcu- 
lated, once the DM profile has been determined with the 
prescription outlined above. In the case of scenario B, 
which leads to higher gamma-ray fluxes than scenario A, 
the differential flux from a mini-spike at distance D can 
be written as 38] 



*{E t D) = *og( 



av 



D 
kpc 

pc 



10 -26 
2 



cm 3 s 1 
P(7s P ) 



(- 



VlOOGeV/ 

2 



10 2 GeVcm- 3 

^"cut 



10" 3 pc 



(6) 



where the radius of the spike is r sp « 0.2^ |74j ]. and 7 S 



with $0 = 9 x 10 _10 cm _2 s _1 . The gamma-ray spectrum 
per annihilation dN/dE depends on the nature of the 
DM particle, and both numerical calculations and ana- 
lytic fits are available in the literature for all possible 
annihilation channels. This formula is valid for mini- 
spikes forming adiabatically from an initial NFW profile, 
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and under the (very good) approximation r sp 3> r cu t- In 
this case, one can easily verify that the "luminosity" of 
a mini-spike in terms of DM annihilation is of the same 
order of magnitude as that of the entire Milky Way halo. 

Although one would naively expect that the fluxes 
scale with av/m x , in the mini-spike scenario the DM 
profile itself depends on m x and av, and the final lu- 
minosity of the objects, for an initially NFW profile, is 



proportional to ~ (o"i>) 



2/7 m "9/7 



E 

o 



uT 

A 

O 

o 

o 



-13 



] 1 1 — 1 1 1 1 1 1 1 1 1 — 1 — 1 


1 1 1 1 1 


1 1 — r-i 


I 1 1 1 IJ 




<o v> 


= 3x10" 26 cm 


v : 






E th = 100GeV - 






E th = 10GeV 








E th = 1GeV 


























— 








] 


; ; ; i 
















-■i 












i i 









10" 



10" 



1 10 
m x [TeV] 



Figure 3: Integrated flux $ 7 expressed in cm" s _1 above the 
energy threshold as a function of the neutralino mass m x for 
thresholds of 1, 10 and 100 GeV, respectively, and a velocity- 
weighted annihilation cross section av = 3 x 10 -26 cm 3 s _1 . 
The quoted error bars correspond to the r.m.s. of the inte- 
grated flux distribution. The maximum flux is obtained for 
neutralino masses well above the energy threshold of the in- 
strument (see text for details). 

Fig. [3 shows the integrated luminosity function of 
IMBHs in a Milky Way-sized halo, i.e. the number of 
black holes Vbh producing an integrated gamma-ray flux 
higher than a given flux, as a function of the integrated 
flux. This can be understood as the number of black 
holes that can be detected with a telescope of given inte- 
grated flux sensitivity. The point-source sensitivity (5a, 
25 hours at 20° zenith) for H.E.S.S. is plotted for com- 
parison. The integrated luminosity is shown in the case 
of the aforementioned scenarios A and B, for three differ- 
ent dark matter particle masses and an annihilation cross 
section av = 3 x 10 -26 cm 3 s _1 . This value allows for the 
thermal relic density of the DM particle to account for 
the measured cold dark matter density Ocdm/i 2 —0.1 
(see Eq. (3.4) of Ref. [|). Here, DM particles are as- 
sumed to be neutralinos annihilating in the bb channel. 
For the assumed value of av, a large number of IMBHs 
in the Milky Way is within the reach of H.E.S.S. for sce- 
nario B. In what follows, we will concentrate on that sce- 



nario. Fig. [3] presents the mean integrated gamma-ray 
flux per IMBH for various energy thresholds as a func- 
tion of the mass of the DM particle annihilating into bb. 
Also displayed are the error bars corresponding to the 
r.m.s. variation of the integrated flux distribution. For 
a 1 GeV threshold, well suited for gamma-ray satellite 
experiments, the maximum flux is obtained for a DM 
particle mass of ~ 80 GeV. This maximum comes from a 
balance between the factor m x 9 ^ 7 and the integral of the 
annihilation spectrum up to the DM particle mass (see 
Eq.[5]). Adopting an energy threshold of 100 GeV, as ap- 
propriate for Cherenkov telescopes such as H.E.S.S., the 
largest fluxes are obtained for a mass of ~ 5 TeV. For 
this mass, the mean of the integrated flux distribution is 
4.5 x 10 -11 cm~ 2 s _1 . For masses close to the experimen- 
tal threshold, the integrated flux increases with the dark 
matter mass. Well above the threshold, the standard 
regime is recovered, with fluxes decreasing with m x . 



III. H.E.S.S. 

A. The H.E.S.S. instrument 

The H.E.S.S. (High Energy Stereoscopic System) ar- 
ray is dedicated to VHE gamma-ray astronomy. The 
instrument is composed of four Imaging Atmospheric 
Cherenkov Telescopes located in the Khomas Highland 
of Namibia at an altitude of 1800 m above sea level. 
This southern location is well suited for observations to- 
wards the inner region of the Galactic halo. Each tele- 
scope consists of an optical reflector of 107 m 2 composed 
of 382 round mirrors [76[. The Cherenkov light emit- 
ted by charged particles in the electromagnetic shower 
initiated by the primary gamma-ray is focused onto a 
camera equipped with 960 photomultiplier tubes (PMTs) 
of 0.16° individual field of view [77|. The total field of 
view of H.E.S.S. is 5° in diameter. The stereoscopic tech- 
nique allows for accurate reconstruction of the direction 
and the energy of the primary gamma-ray [78|. H.E.S.S. 
has an angular resolution of about 5 arc-minutes and 
a source location accuracy of ~ 30" for strong sources. 
The sensitivity for point-like sources reaches 2 x 10~ 13 
cm - ^" 1 above 1 TeV for a 5a detection in 25 hours at 
20° zenith [z|. 



B. Observations and data analysis 

The data used in this analysis were collected between 
2004 and 2007 during the Galactic plane survey with the 
four telescope H.E.S.S. array. The Galactic plane has 
been observed between ± 3° in latitude and from -30° 
to 60° in longitude relative to the Galactic Center in 884 
pointings. Runs are taken mainly with 28 minute dura- 
tion at pointing positions with a typical spacing of 0.4° in 
longitude and 1° in latitude. For the analysis described in 
this paper, the observations dedicated to known sources 
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at other wavelengths are excluded. Astrophysical models relation: 
of TeV emission are available for these sources so that 
the search is focused on regions where no standard astro- 
physical emitters have been detected by H.E.S.S. After 
the standard quality selection procedure [79j and dead 
time correction, the data set amounts to ~ 400 hours of 
live time and a mean zenith angle of all observations is ~ 
30° resulting in a typical energy threshold of 200 GeV. 

Following the standard calibration of the shower im- 
ages from PMT signals [SO ], the event reconstruction 
scheme using the combined Model-Hillas analysis is ap- 
plied to the data to select the gamma events. The Hillas 
reconstruction is based on the Hillas geometrical mo- 
ment of the image 81]. The Model analysis is based 
on a pixel-by-pixel comparison of the image to a tem- 
plate image generated by a semi-analytical model of the 
shower [82| . Both methods yield a typical energy resolu- 
tion of 15% and an angular resolution at 68% full contain- 
ment radius better than 0.1°. The data analysis is done 
with a combination of these two methods to improve the 
hadronic background rejection (8l| . An additional cut 
on the primary interaction depth is also used to improve 
background rejection. After a cleaning of the images, the 
direction, energy, impact parameter and primary inter- 
action point are reconstructed for each gamma event. A 
cut on the image size at 60 photoelectrons is used to ob- 
tain better sensitivity to weak gamma-ray excesses. The 
background level is estimated using the template model 
method [83[ as described below. This allows to estimate 
the background level at each sky position. 



C. Gamma-ray sensitivity map 

In order to study the H.E.S.S. sensitivity in the Galac- 
tic plane survey to dark matter annihilations from mini- 
spikes, we perform the analysis on the Galactic plane sur- 
vey data. The H.E.S.S. flux sensitivity map is derived us- 
ing the following sky maps in Galactic coordinates which 
are estimated within the survey range: 

• the observed gamma-ray map, 

• the normalized measured background map, 

• the effective time exposure map (in seconds), 

• the gamma-ray acceptance map (in m 2 ). 

The gamma-ray map is obtained after the event selection 
and reconstruction procedures described above. The sky 
map is divided into bins of 0.02° x 0.02° size and an 
oversampling radius of 0.1° is applied to the data. The 
oversampling smoothes the map with a top-hat function 
of radius of 0.1° to match the point spread function of 
the instrument. In each bin of the map in Galactic co- 
ordinates, the background rate is estimated by the tem- 
plate background method (see [83| for details) and the 
normalized background map is obtained by the following 



(b,l) = N Bkg (b,l) x 



^4cc 7 (6, 1) 
Acch(b, I) 



(7) 



where the acceptances for gamma-ray- like events, Acc^, 
and hadron-like events, j4cch, are computed in each po- 
sition of Galactic longitude b and latitude I taking into 
account collection area corrections. Using the observed 
gamma-ray map and the normalized background map, we 
derive a sensitivity map expressed in terms of an upper 
limit (at 90% C.L) on the number of gamma events above 
the nominal energy threshold in each sky position. Using 
the exposure and acceptance maps, the flux sensitivity 
map for gamma-rays from dark matter annihilation is 
then determined by : 



90%C.L 
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(M) 



N*> %aL (b,l) 



DM dN 

Lth dE 1 



dE~, 



T obs JO 



A eS {E 1 ,z{b,l)A^l))^{E 1 ) 



.(8) 



dE-y dr 



Here, A e g is the effective area for gamma-rays which is 
a function of the gamma energy E 7 and the sky position 
(b,l). The dr integration denotes the averaging accord- 
ing to the zenith angle z and offset 9 distribution of the 
observation livetime. The annihilation spectrum of the 
dark matter particle of mass todm is denoted by dN/ dE 7 
and is integrated from the nominal energy threshold E t h, 
which is about 100 GeV in the survey, up to mdm- The 
integral over r is calculated over the total observation 
livetime T ^ s . 

Fig. 0] shows the experimentally observed sensitiv- 
ity map in the Galactic plane from Galactic longitudes 
I = -30° to I = +60° and Galactic latitudes b = -3° to 
b = +3°, for a dark matter particle of 500 GeV mass an- 
nihilating into the bb channel. The H.E.S.S. sensitivity 
depends strongly on the exposure time and acceptance 
maps which are related to the choice of the pointing po- 
sitions. The flux sensitivity varies along the latitude and 
longitude due to inhomogeneous coverage of the Galactic 
plane. In principle, the sensitivity map depends on the 
dark matter annihilation spectrum. However, as can be 
seen from Eq. ([5]), the spectrum is balanced by the effec- 
tive area which mainly drives the result of the integral. 
The particle mass is not expected to bring about strong 
variations in the map as long as the mass is larger than 
100 GeV. 

In the band between —2° and 2° in Galactic latitude, 
a DM annihilation flux sensitivity at the level of 10~ 12 
cm _2 s _1 is achieved. H.E.S.S. thus reaches the required 
sensitivity to be able to test dark matter annihilations 
from mini-spikes in the context of one relatively favorable 
scenario for IMBH formation and adiabatic growth of the 
DM halo around the black hole. 

Deeper observations of the GC and at Galactic longi- 
tude of ~ —20° allow the flux sensitivity of ~ 5 x 10~ 13 
cm _2 s _1 for a 500 GeV DM particle annihilating in the 
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Galactic longitude [deg.] 

Figure 4: H.E.S.S. sensitivity map in Galactic coordinates, i.e. 90% C.L. limit on the integrated gamma-ray flux above 100 
GeV, for dark matter annihilation assuming a dark matter particle of mass m x = 500 GeV and annihilation into the bb channel. 
The flux sensitivity is correlated to the exposure and acceptance maps. In the Galactic latitude band between -2° and 2° , the 
gamma-ray flux sensitivity reaches 10 -12 cm -2 s -1 . 



bb channel. For 6 = 0° and I = —6°, the flux sensitivity is 
~ 1CP 12 cm _2 s . For |6| > 2°, the sensitivity is deterio- 
rated due to a weaker effective exposure. For 6 = 0° and 
/ = —0.5°, near the Galactic Center, the flux sensitiv- 
ity is ~ 10~ 13 cm~ 2 s -1 in the 100% BR bb annihilation 
channel and ~5x 10 -14 cm _2 s _1 in the t + t~ channel. 



IV. RESULTS AND DISCUSSION 

H.E.S.S. observations (2004-2006) of the Galactic 
plane allowed to discover more than 20 VHE sources • 
Some of them have been identified owing to their coun- 
terparts at other wavelengths, but almost half of the 
sources have no obvious counterpart and are still uniden- 
tified [85]. For these sources, an accurate reconstruction 
of their energy spectra has been carried out. All spectra 
were consistent with a pure power-law of spectral indices 
between 2.0 and 2.5, spanning up to two orders of mag- 
nitude in energy above the energy threshold, as shown in 
Fig. 8 of Ref. [85| . None of them exhibits an energy cut- 
off, characteristic of dark matter annihilation spectra, in 
the energy range from ~ 100 GeV up to 10 TeV. Further- 
more, the detailed study of their morphology 85] show 
that all the sources have an intrinsic spatial extension 
greater than ~5', while mini-spikes are expected to be 
point-like sources, since the bulk of the gamma-ray emis- 
sion comes from a region of size « 10~ 3 pc at a typical 
distance from GC of 10 kpc. In the survey region dis- 
cussed here, only three point-like gamma-ray sources are 



detected: HESS J1826-148, HESS J1747-281 and HESS 
J1745-290. They are identified with well-known objects: 
LS 5039, GO. 9+0.1 and the Galactic Center, respectively. 
H.E.S.S. has detected so far no IMBH candidate within 
the survey range. 

Based on the absence of plausible IMBH candidates in 
the H.E.S.S. data, we can derive constraints on one sce- 
nario for neutralino or LKP dark matter annihilations. 
These constraints are shown as upper limits on the anni- 
hilation cross section (Figures [5] and [7]) but are actually 
constraints on the entire gamma-ray production scenario, 
assuming the mass function of scenario B and a coreless 
halo model. In the survey range, the expected mean num- 
ber of IMBHs is 4.3, within the assumptions of scenario 
B, as shown in Fig. [5] The number of IMBHs in the whole 
halo depends on the assumed cosmological parameters in 
the formation scenario [86], in particular the reioniza- 
tion rcdshift. The value of 4.3 unmerged IMBHs is for 
a reionization redshift of 16. The redshift of 10.8 ± 1.4 
determined from latest WMAP data [87[ favors a larger 
number of IMBHs since they had less time for merger 
events to occur [86] . The distribution of the IMBH num- 
ber has an intrinsic rms of 2.3. The probability to find 
no IMBH in the survey range is 2.5%. The simulations 
of the Milky Way halo show that the IMBH number dis- 
tribution in the sky survey range is well fitted with a 
Poisson distribution with mean 4.4. 

For a given gamma-ray flux sensitivity <E>, we define 
the probability density function of black holes yielding 
an integrated gamma-ray flux larger than <E>, cL 2 Nbh{> 
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masses in the TeV energy range, we obtain limits on one 
mini-spike scenario (scenario B). The limits on av are at 
the level of 10 -28 cm^ 3 s~ 1 for the bb channel. A rapid 
decrease in sensitivity is observed for WIMP masses less 
than ~ 800 GeV. This corresponds to the threshold effect 
seen in Fig. [31 Predictions for SUSY models are gener- 
ated using the DarkSUSY code [88| in a phenomenologi- 
cal MSSM (pMSSM) framework characterized by the fol- 
lowing independent parameters: the common scalar mass 
mo, the gaugino mass parameter M2, the higgsino mass 
parameter /1, the CP-odd Higgs mass Ma, the trilincar 
couplings A t and the Higgs vacuum expectation value 
ratio tan (3. The set of parameters for a given model is 
randomly chosen in the pMSSM parameter space encom- 
passing a large class of pMSSM models, as described in 
Table [IJ Models providing a neutralino thermal relic den- 
sity, fiDM^ 2 , in the range [0.08,0.12] are overlaid to ac- 
count for the cold dark matter density inferred from the 



Figure 5: Distribution of the number of IMBHs for the sce- 
nario B in the survey range corresponding to a field of view 
between ±3° in Galactic latitude and from -30° to 60° in 
Galactic longitude. The distribution is obtained from 200 
stochastic realizations of Milky Way-like halos (see Sec. II. A 
for details). The mean of the distribution is 4.3 with an rms 
of 2.3. The distribution is well fitted to a Poisson distribution 
with mean 4.4 (black dashed line). 

$)/dbdl, from the IMBH catalog used in [38] by : 
N B a (detected, > $) = 

f b=+3° f l=+60° d 2 N (> $ n 
Jb=-3° J2 = -30 o aoal 

Nbh (detected, > $) denotes the expected number of 
black holes yielding a gamma-ray flux larger than $ in 
the survey, $ corresponding to the H.E.S.S. sensitiv- 
ity. The total number of black holes in the survey is 
computed by integrating the probability density function 
d 2 NBn{> <&)/dbdl over the latitude and longitude ranges 
of the Galactic plane survey. Since no H.E.S.S. source is 
a plausible IMBH candidate, we thus calculate for each 
dark matter particle mass the limit at 90% C.L. on the 
velocity- weighted cross section av, assuming a Poisson 
distribution. 

Fig. O shows the exclusion limit at the 90% C.L. on 
av as a function of the neutralino mass todm ■ The neu- 
tralino is assumed to annihilate into bb and t + t~ with 
100% BR, respectively. Below ~ 2 TeV, the upper con- 
tour of the gray shaded area is given by the bb assumption 
and yields the conservative exclusion limit. Above 2 TeV, 
the upper contour corresponds to the t + t~ annihilation 
spectrum. The maximum sensitivity of H.E.S.S. to DM 
annihilation for the bb channel is achieved for masses of 
order 5 TeV as expected from Fig. [3J For neutralino 




Figure 6: Constraints on the IMBH gamma-ray production 
scenario for different neutralino parameters, shown as upper 
limits on the annihilation cross section av as a function of the 
mass of the dark matter particle wdm, but with a number 
of implicit assumptions about the IMBH initial mass func- 
tion and halo profile (see text for details). For scenario B, 
the probability of having no observable halos in our Galaxy 
is 10% from Poisson statistics making these limits essentially 
90% confidence level exclusion limits for this one particular 
(albeit optimistic) scenario (grey shaded area). The dark mat- 
ter particle is assumed to be a neutralino annihilating into bb 
pairs or r + r~ pairs to encompass the softest and hardest an- 
nihilation spectra. The limit is derived from the H.E.S.S. flux 
sensitivity in the Galactic plane survey within the mini-spike 
scenario. SUSY models from the pMSSM (black points) are 
plotted together with those satisfying the WMAP constraints 
on the dark matter particle relic density (magenta points). 
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Table I: pMSSM parameter space randomly scanned to gen- 
erate SUSY models. A set of parameters corresponds to a 
specific pMSSM model. 



Parameter Minimum Maximum 



m 1 TeV 30 TeV 

M 2 1 TeV 50 TeV 

H 1 TeV 50 TeV 

m A 1 TeV 50 TeV 

A t , b -300 GeV +300 GeV 

tan/3 3 60 



channel. The predictions for av are parameterized using 
the results of [89j where it is derived in the case of the lat- 
ter assumption on the annihilation scheme. The AMSB 
models fulfilling in addition the WMAP constraints are 
overlaid. Models yielding neutralino masses between 900 
GeV and 6 TeV are excluded including those having 
a neutralino thermal relic density compatible with the 
WMAP measurements. Predictions on the annihilation 
cross section for pairs of from UED theories [90] are 
also plotted as well as those satisfying the WMAP con- 
straints. In this scenario, masses in the range from 
0.8 to 6 TeV can be excluded. 



measurements of the CMB anisotropies of the Wilkinson 
Microwave Anisotropy Probe (WMAP) satellite. Some 
models in the mass range 0.8 - 6 TeV can be excluded. 
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Figure 7: 90% C.L. exclusion limits on the annihilation cross 
section av as a function of the mass of the dark matter par- 
ticle rnoM for a Kaluza-Klein boson B"' (solid blue line) and 
a pure wino W in AMSB (solid red line). The limit is derived 
from the H.E.S.S. flux sensitivity in the Galactic plane survey 
within the mini-spike scenario. Kaluza-Klein models in UED 
scenario (blue dashed line) and AMSB models (red dotted 
line) are plotted together with those satisfying the WMAP 
constraints on the dark matter particle relic density (pink 
dashed line and green dotted line, respectively). 

Fig. [7] shows the 90% C.L. exclusion limit on av as 
a function of the DM particle mass todm in the case of 
AMSB and Kaluza-Klein scenarios. In both cases, the 
BRs of each annihilation channel entering in the calcu- 
lation of the total annihilation differential spectra are 
uniquely determined as discussed in section II. In the 
AMSB scenario, the neutralino is considered here to be 
a pure wmo annihilating with 100% BR into the W + W~ 



V. SUMMARY 

Observational clues for the existence of IMBHs start 
to accumulate. If they indeed exist, IMBH could be dark 
matter annihilation boosters. The prospects for detect- 
ing dark matter annihilation around IMBHs have been 
widely discussed in the literature [H, [H, H^, [9l|, [51, HU . 
In this work, we derive the first experimental constraints 
on the mini-spike scenario of Ref [38| . 

Using H.E.S.S. data collected in the Galactic plane 
survey, we show that H.E.S.S. has the required sensitiv- 
ity to probe gamma-rays from dark matter annihilation 
in mini-spikes around IMBHs believed to populate the 
Milky Way halo. The new analysis using ~ 400 hours of 
data taken in the Galactic plane but not foreseen initially 
for this purpose, allows to derive flux sensitivity limits for 
indirect dark matter search. Combining all the survey 
data, the gamma-ray flux sensitivity map is derived for 
dark matter annihilation in the region [—30°, 60°] in lon- 
gitude and [-3°,3°] in latitude. We show that H.E.S.S. 
reaches a flux sensitivity of ~ 10~ 12 cm~ 2 s above 100 
GeV. 

For the first time, a dumpiness scenario has been 
tested in a large field of view with an IACT. Strong 
constraints are obtained in one of the two IMBH for- 
mation scenarios discussed 3S\ (scenario B). The ab- 
sence of plausible candidates for Galactic IMBHs in the 
H.E.S.S. Galactic plane data set allows us to put con- 
straints on one of the more optimistic scenarios for de- 
tecting neutralino or LKP annihilation from mini-spikes 
around IMBHs. The first experimental exclusion limits 
at 90% C.L. on the velocity-weighted annihilation cross 
section as a function of the dark matter particle mass 
within the mini-spike scenario are derived. Predictions 
from various WIMP particle physics scenarios are con- 
strained. 

Since the characteristic annihilation flux from a IMBH 
varies as (av) 2 / 7 , reflecting the depletion of the inner part 
of dark matter halo due to annihilation during the life- 
time of the IMBH, limits derived on av in the absence 
of a detection are proportional to (^min/^rof) 7 ^ 2 , where 
^min is the sensitivity limit and <j> re f the predicted an- 
nihilation flux for a nominal value of av. Uncertainties 
in $ re f, reflecting for example the imperfect knowledge 
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of the exact shape of the dark matter spike close to the 
IMBH, cause correspondingly enlarged uncertainties on 
the limits on av. Similarly, variations in the predicted 
number and mass of IMBHs may influence the limits; for 
example, at least about 50 IMBHs have to be contained 
in the Galactic halo in order to set a meaning flux up- 
per limit, given the limited solid angle coverage of the 
H.E.S.S. survey. The limits given here apply within the 
formalism and approximations of scenario B of [38| . 

The analysis described here can be adapted to other 
particle physics parameters, e.g. different types of DM 
particles and annihilation channels. Constraints are 
derived within the mini-spike scenario although the 
method developed here is generic and is suited whatever 
the assumed dark matter clump scheme. The flux 
sensitivity map is a powerful tool to investigate the 
sensitivity to other t ypes of overdensities such as small 
scale clumps [13, Ha . I94L 1951 o r compact dark matter 
structures like spikes [25l l75l. 1 9(1 l97l . [§8| . 
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